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ABSTRACT 



Averaged spin-spin correlation function squared < o{Q)o{R) is calcu- 
lated for the ferromagnetic random bond Potts model. The technique being 
used is the renormalization group plus conformal field theory. The results are 
of the e - expansion type fixed point calculation, e being the deviation of the 
central charge (or the number of components) of the Potts model from the 
Ising model value. Calculations are done both for the replica symmetric and 
the replica symmetry broken fixed points. The results obtained allow for the 
numerical simulation tests to decide between the two different criticalities of 
the random bond Potts model. 
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1 Introduction. 



For the ferromagnetic random bond Potts model there exist two different fixed point 
solutions, one which is replica symmetric [1,2] and a second in which the replica symmetry 
is broken [3] . Both fixed points are of a universal nature and one or another is reached 
depending on the initial conditions for the coupling constant in the renormalization group 
(RG) equation. More specifically, if one starts with the model defined on a lattice, then 
the partition function shall be given by: 

Z(/3)=^exp{-/3//[a]} (1.1) 
{^} 

where H[a] is a nearest neighbor interaction hamiltonian for classical spins {ax}: 

H[a]^J2J-,o.V{ax,ax+a) (1-2) 

x,a 

Each spin is taking q values, in case of g-component Potts model; x stands for lattice 
sites, {x,a) - for lattice bonds: in case of a square lattice a — 1,2; V{ax,crx+a) is the 
spin-spin interaction potential: in case of the Potts model the usual choice would be 

V{a,(7') oc 1 — 6ay. {Jx,a} sxB the bond coupling constants. They are supposed to be 
ferromagnetic but taking random values, independently for each lattice bond, with some 
distribution, characterized by a width go, around a mean value Jq. 

Weak disorder corresponds to a small value of go : go ^ Jo- In this case the model 
could be studied in the continuum limit, if (3 is close to a critical value (3c, {/3 — f3c)/(3c = 
T <^ 1. In this limit the partition function could be given in the form: 

Z{/3) = Trexp{-Ho-H,} (1.3) 

Tr stands symbolically to represent the sum over the spin configurations, but in the 
context of the continuum limit theory. Its explicit realization is not important because 
we shall eventually be dealing with correlation functions, and these are defined unam- 
biguously by the corresponding conformal field theory. Ho stands for a hamiltonian, or a 
field theory action, of the conformal field theory corresponding to a given g-states Potts 
model, being defined on a perfect lattice with the spin-spin coupling constant Jq for all 
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lattice bonds and taken at its critical point, (3 — (3c] Hi represents a deviation from the 
critical point and it contains disorder. Hi could be given in the following form: 

Hi^ J (fxr(x)e(x) (1.4) 

where 

t{x) oc pj{x) - PcJo (1-5) 

is the random temperature parameter of the continuum limit theory; x takes values on the 
continuum 2D plane; J^^a the lattice is replaced by J{x); e{x) is the energy operator 
of the Potts model replacing V{ax,(Tx+a) on the lattice. As far as critical properties are 
concerned the reduced continuum limit form of the model, defined by the eqs. (1.3), 
(1.4), is sufficient. For simplicity we shall assume that t{x) has a Gaussian distribution, 
for each x, with a width qq, so that 



t(x) = To = (/3 - /3e)//3c (1.6) 

{r{x)-To){T{x')-To) = 50^(2) {x - x') (1.7) 
The partition function (1.3) is of the form: 

Z{(3) = Tr exp{-/7o - j d^xT{x)e{x)} (1.8) 

To take the average over the disorder one introduces replicas, n copies of the same model: 

n „ n 

{Z{P)Y = Tre^vi-Y^H^^ - / d?xT{x)Y.ea{x)} (1.9) 

o=l ■' a=l 

and then one takes the average: 

n „ n „ n 

{Z{(3)r = Trexp{-J2Hl,"^-ro d^xj^ Saix) + go d\Y.^a{x)e,{x)} (1.10) 

a=l ^ a=l ^ aj^b 

One arrives in this way at a homogeneous field theory of n coupled models with the 
coupling action: 

Hi^i^-go [d'x^e^:h,{x) (1.11) 
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In the non-diagonal terms only, a ^h, are being kept. The diagonal ones could be 
put back into Z]a=i -^o"'*- Moreover, in case of the Potts model the energy operator e{x) 
corresponds to the operator $1^2 (a^) of the corresponding conformal theory. Its operator 
algebra is known 

^1,2^1,2^^1,1 + ^1,3 (1-12) 

Here $1^1 = J is the identity operator, and the operator $1 3(3;) is irrelevant, 2Ai 3 > 2. 
So there will be no relevant diagonal subtractions from Z^ab^a^b, the diagonal terms 
can just be dropped. In our analysis of the random bond Potts model we shall take 
the replicated field theory form of it, equations (1.10), (1.11), as our starting point. 
In the renormalization group analysis, with the interaction term given by the equation 
(1.11), the qualitatively different initial conditions for the renormalization group equation 
correspond to cither taking the initial coupling constant as it is in eq.(l.ll), g{!^ = ^) = go 
(^ is the RG parameter), or to breaking the replica symmetry initially by putting H-^^ii 
into the form: 

^int = ~ / 'f'x^gahea{.x)ei,{x) (1.13) 
This corresponds to assuming different couplings, initially, for different replicas. 
Replica symmetric form of gah is 

9ah = 5-0, all a^h (1.14) 

and gah = for a = 6. Taking g^h with different components corresponds to replica 
symmetry breaking perturbation. Relevance of this type of perturbation has been ob- 
served initially in [4] for the random XY model and in [5] for the standard 99^ theory, 
with randomness. These two models move, under RG, into a strong coupling regime. 
The first fixed point with a replica symmetry breaking has been found in [3] for the 2D 
random-bond Potts model. For the moment this remains to be the only solution of this 
type. 

To repeat it again, the random bond Potts model, i.e. the model with iJj^^ in (1.13), 
has two different fixed point solutions: 

1) replica symmetric, if initially gah of the form in eq.(1.14); 
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2) replica symmetry broken, if initially Qab is different from (1.14). 

To check which one of the two solutions realizes for the original spin model defined on 
a lattice with random bonds, eqs.(l.l), (1.2), in the numerical simulation experiment in 
particular if not in the real one, we shall calculate in this paper one particular observable 
quantity for which the replica symmetry breaking effects are expected to be most pro- 
nounced. The calculations will be done for both RG fixed points to make the verification 
problem well defined. 

The quantity in question is the averaged spin-spin correlation function squared 

< a(0)a(R) >2 (1.15) 
The related quantity is the renormalization group amplitude 

for the replicated model operator 

Oabix) = aa{x)ab{x) (1.17) 

Here the spins Ua, ab, belong to different replicas, a b, and are taken at the same 
point. The amplitude (1.16) is the spin-spin "overlap function" for our critical model, if 
one uses the terminology of the theory of spin glasses [6] . 

Further discussions will be postponed until the formal analytic solutions are obtained. 

In the next section the analytic problem will be defined and the RG equation for the 
amplitude .^06 (0 will be obtained, up to second order in the coupling constant QabiO- 
Section 3 the amphtude Zab and the correlation function < a{0)a{R) >^ will be calcu- 
lated. We shall also define the corresponding magnetization. In Section 4 the results will 
be analyzed, in particular as applied to the 3 and 4 component Potts models. We shall 
also consider in this section the distribution function for products of local magnetizations. 
This distribution is encoded in the spin-spin "overlap function", the amplitude Zab, in a 
way analogous to that for the corresponding object in the spin-glass theory. In section 
5 results of numerical simulations will be presented for the 3 and 4 - component Potts 
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models. Section 6 is devoted to conclusions and discussions. In the Appendices the cal- 
culations of particular complicated integrals are exposed, which appear in the calculation 
of coefficients of the RG equation for ZabU)- 

2 Definition of the problem and calculation of the 
RG equation for Zab- 

The RG equation for the coupling constant Qab in the action H^^^, eq.(1.13), has been 
derived and solved in [3], for the replica symmetry broken fixed point. The replica 
symmetric solution has originally been found in [1] . We shall use those results but first 
we shall concentrate on the RG equation for the amplitude ^a&(0- This function occurs 
in the calculation of the correlation function < a{0)a{R) >^. 

In terms of replicas the averaged spin-spin correlation function, {a{0)a{R))'^ 
— lim„_^o ^^^zi) Sa7^&(<^a(0)cr(j(it!)(T6(0)(T6(it!)), can be represented in the following form 
(see Sec.3, eqs.(3.7)-(3.9)): 

< a{0)a{R) >^ ^ lim J_ < ^ a„(0)a5(0) ^ <7e(it:)<7,(it:) > (2.1) 
Then the operator to be renormalized is: 

Oab{x) = aa{x)ab{x), aj^b (2.2) 

As we are going to use the perturbative RG, we have to consider the exponential of H^Yrt, 
the way it would enter the partition function (1.10), in the presence of the operator 
Oab{x): 

Oab{x) exp{-ifj^^} = aa{x)ab{x) exp{ / (fy ^ gab£a{y)£b{y)} (2.3) 

a^b 

We expand exp{— iJjj^^}: 
Oab{x){l - i/int + ^(^int)' + •••) = Oabix) - Oabix)H^nt + ^0.,(a;)i/f^^ + • • • (2.4) 

and then we do all possible contractions and operator algebra, which lead to reproducing 
the operator Oab{x). Contractions in the case of our non-gaussian field theory amounts 
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to equalizing replica indices in all possible ways which occur under the summation over 
indices and then replacing the products of operators with the same index by correlation 
functions (or using directly the operator algebra) of the unperturbed theory. We remind 
that the unperturbed theory is a collection of independent and identical Potts models. 
Each Potts model, in the continuum limit, is a minimal conformal field theory with 
the energy operator e{x) represented by the primary field $1^2(2^), and with the central 
charge related to the number of components q of the original Potts model as defined on 
the lattice. We shall specify this relation later, see also [7]. 

Proceeding in this way, one gets in the first order: 

-Oab{x)H^^^ = aa{x)ab{x) d^yY,gcd£c{y)£d{y) 

'^9abcra{x)crb{x) / (fy-: (2.5) 

Jl<\y-x\<a \X — y\^^^ 

We have contracted here the indices b and c, and a and d (or b and d, and a and c, which 
amounts to a combinatorial factor of 2), and then we have used the operator algebra: 

a(x)e(y)= ^ a(x) + ■ ■ ■ (2.6) 

\x — y\^^ 

D stands for the operator algebra coefficient D"^^. For minimal conformal theories these 
coefficients have been calculated in [8] . 

Integration in (2.5) is around x, over the distances ranging between the "old" and 
the "new" cut-offs. The old one we choose to be equal to 1 (like it would be on a lattice, 
when distances are measured in lattice spacings) and the new one a, a >> 1. 

e{x) is the primary field $1^2, so one has: 

A, = Ai,2 + Ai,2 = 2Ai,2 (2.7) 

We use next the Kac formula for conformal dimensions of primary fields $„/^„(x): 

^ ^ ^ («n' + a+nf - {a- + a-f 

The parameters «+, «_ are related to the central charge of the Virasoro algebra: 

c = l-24a^. 



a± = aQ±Jl + al (2.9) 



Notice that 

q;+q;_ = -1 (2.10) 

We shall use in the following the parameter a'\. to specify the models, instead of the 
central charge c. For Ising model c — \, oi\ — ^. One gets conformal theories for the 
critical Potts models with the number of components q varying continuously from 2 to 
4 if the central charge is taken to vary between 1/2 and 1, or varying between 4/3 
and 1 [7]. For the Ising model Ai_2 = |, = 1, and the perturbation i^jj^t (1-13) is 
marginal, gab is dimensionless. This perturbation becomes slightly relevant if, following 
Ludwig [1], one takes 

4 = ^-e (2.11) 
and one studies the Potts models by the e-expansion RG assuming e to be small. 
The dimension of the energy operator is now given by: 

A.= 2A.,= '"-+^°^)'-'"-+"^''=l-^e (2.12) 

The first order correction in (2.5) takes the form: 



1 

'2gaj,aa{x)a^{x)D'' [ — ^ = 2gai,<Ja{x)ab{x)D''2n-a^' (2.13) 
Ji<\y\<a yr '^'^ 3e 



n<\y\<a \y\ 

In integrating over scales from 1 to a one gets a factor ^{a^^ — 1). We have replaced it, 
in a standard way, by \a^^, assuming a^^ ^ 1. 

Next, if one introduces the amplitude Zab and studies renormalization of the operator 

Oabix) = Zabcraix)ab{x) (2.14) 

then (2.13) will correspond to the first order correction to Zab- 

Szj^^ = Zab^TrD^gab^-/' (2.15) 

This correction, the corresponding RG equation, and the renormalized amplitude Zab has 
first been defined in [1], (for the replica symmetric case), as well as the corresponding 
amplitudes for higher moments of a: Oab-d = 'Z'ab-d'^a^b • • • Cd, a 7^ 6 7^ ■ ■ ■ 7^ d. 
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It has been shown in [3] that the rephca symmetry breaking effects appear in the 
second order. One needs at least two orders, for all the quantities, to treat these effects. 
For this reason our analytic problem will be to define renormalization of Za\, in two orders. 

Prom eqs.(2.3), (2.4) one gets in the second order: 

Oah^iiiYafd^ = Zab(Ta{.x)ab{x)]- j d^yY,9cd£c{y)£d{y) f d:^y'Y,gef£e{y')£f{y') 

^ Df^ + + (2.16) 

We have defined as Df^ the diagrams of the second order. The first diagram is of the 
form: 

-Of^ = ^Zab(Ta{^)(Tb{^)\Y^9bdgad 

^ d 

X jd^yj < a{0)e{y)a{oo) >< a{0)e{y')a{oo) >< e{y)e{y') > 

= 4Z„,a„(0)<7,(0)(^^)„, / A / ^V^ i^^^^ \y-y'r^ ^^"^^^ 

To simplify somewhat the expressions in the integrals we have put the operator Oah{x) = 
Zabcra{x)ab{x) at a; = 0. 

To get D^^ we have made equal in (2.16) b — c and a — e (we remind that a ^ h, 
c ^ d, e ^ f : diagonal elements of Zab,gcd, Qef are assumed to be zero). Next we 
have used the operator product expansion (2.6) for the products a{{))£{y) and a(0)£(y'), 
keeping just the first term, the only relevant one, or, which is the same, we have calculated 
the projections of (T{Q)e{y) and of (T{Q)e{y') on the spin operator placed at infinity. In 
addition, we have assumed the replica coupling matrix Qcd to be symmetric, Qd — ddc, as 
this is the case for the Parisi matrices and for the fixed point solution for gab in [3]. In 
this case J^ddbddda = Y^dQadQdb = {g^)ab- Piually, the extra factor of 8 in the first line in 
(2.17) is due to combinatorics: there are 8 ways to make coupling of indices which give 
the diagram 

The calculation of the integral in (2.17) is straightforward, it is exposed in the Ap- 
pendix A.l. One gets the following result: 

^ = Zabaa{0)abiO)87r'D\g^)abil + g^^'' + ^(1) (2.18) 

9 



where K — 6 log 2. In (2.18) we have kept only the terms which are singular in e, ~ ^ 
and -, the ones which are relevant for the RG. 

The next diagram is obtained from (2.16) by setting h — c — e and, separately, d — f, 
plus all equivalents. This gives: 

oc Z,,aa{0)abmg\b J d^V j < a{Q)e{y)s{y')a{:^) >< e{y)e{y') > (2.19) 

This diagram does not produce singularities in e, it is finite for e ^ 0, see Appendix A. 2. 
As a result, it does not contribute to the RG evolution of Zab- 

Next diagram, D\ \ is obtained by setting h — c — e and a — d — f, plus all 
equivalents. This gives the following expression: 

= iZabaa{0)ab{0)^{gabf J d% J d%'{< a{0)e{y)e{y')a{oo) >f (2.20) 

Calculation of this integral, which is fairly complicated, is described in the Appendices 
A and B. Finally one gets the following result: 

D^s^ = ZaWa{0)abm9ab?{S7T'D'-^^ - (2-21) 

Putting (2.18) and (2.21), for and D^^\ into (2.16), and dropping eq.(2.19), 
one derives the second order correction to the amplitude Zab- 

6Z^^ = Zab{8n'D'{l + eK){g\,^ + Sn'O^g^bf^ - y (2-22) 

Together with 5Z^J in (2.15) one obtains, up to second order: 

Zab = Za, + 5Z^^+SZ^^ 

= Zab{l + ATTD^gabj-^a^' + 87rD2(l + eK){g'')ab-^a^' 

+8n'D\gabr^,a'' - y (^af)'^«'^} (2-23) 

Further in this Section we shall denote by tilde, like Zab or gab, the corresponding renor- 
malized quantities. 
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One obtains, in a standard way, the RG equation for Zab by taking a derivative with 
respect to ^ = log a: 

dZnh dZa 



Jab u.^ab 
= a 

da 

= Z,,{AT,D^gg,a^' + %t,^D\1 + eK){g''U^a'' 

+%^D\g^,f^a'^ - ^^{gabfo^^} (2.24) 

In the r.h.s. of this equation we have to replace Zgjy and gah by the corresponding quanti- 
ties renormalized in the first order. One easily checks, by using the technique described 
above, and, in addition, in case of gabi a dilatation of coordinates, in order to return to 
the cut-off scale a = 1, that to the first order: 

~9ab = 0?'{gab + 47r((7')ab|^) (2.25) 
Zab = Zab{l + ^T^D^gab'^) (2.26) 

We have not been adding the dilatation term in case of renormalization of Zab renor- 
malization because, by its usual definition, we have put in this amplitude only the terms 
which are produced by interactions. The trivial scaling factor in the renormalization of 
the operator Oa&(^) will be supplied later when we shall be calculating the correlation 
function. 

Inversely, in the first order: 

gab = a-''{gab - 47rD\~g\b^J (2.27) 

Zab = Zab{l - ^^D^~9ab^^ (2-28) 

Substituting these expressions in the r.h.s. of (2.24), and keeping terms up to second 
order in g^^, one obtains: 

^ = Za,{AnD'~gab + ^-^D'K{g\b - "^{ga,?} (2.29) 
dt, 6 6 

This is the RG equation for the amplitude Z^b that we were aiming at. 
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It is easy to check that in the rephca symmetric case: 



Qab ^ 9 (2.30) 

Zab(^a(^b Zaa(Jb (2-31) 

- for all a,h {a^h), the RG equation (2.29) takes the form: 

^ = Z{^^D^~g + ^D^Kfin - 2) - (2-32) 

Here n is the number of replicas, to be put equal to eventually. 

In the next Section we shall derive solutions for both equations, (2.29) and (2.32), in 
order to have predictions for both fixed points, symmetric and non-symmetric one. 



3 Solution of the RG equation for the amphtude Zab- 



Correlation function < a{G)a{R) 



We shall drop tildes in c/ab and Zab in the following, as all the quantities that we shall 
use will be the renormalized ones. Also, to simplify somewhat the equations, we shall 
change the normalization of gab' 

- £ (3.1) 

And we shall substitute the value of the operator algebra constant: 

D = Die) ^ Dl^^ie) = \ + 0{^) (3.2) 

It is well known that for the Ising model D^ ^ = 1/2. One could check that the correction 
is ~ , by using the formulas derived in [8] . For our present calculation the correction 
is irrelevant, so we can put the Ising model value D — 1/2. 

The equation for Zab (2.29) takes the following form: 

dZab{i) 



^a6(07a6(0 (3.3) 



lab{i) = \9ab{i) + Y2^9\i))ab - ^^{9ab{i)f (3.4) 
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We remind that the constant K — 6 log 2. For the rephca symmetric case we shall have: 



miii) (3.5) 



7(0 = \m + ^9\i){n - 2) - 1/(0 (3.6) 



We are interested in the correlation function < cr(0)cr(i?) >^ which is expressed through 
replicas by eq.(2.1). Assuming the RG evolution from the lattice cut-off up to the scale 
~ R, one gets: 

X < a„(0)a,(0)ae(l)ad(l) > (3.7) 

Here — log i?; A^°) is the unperturbed dimension of the spin operator in the Potts 
model. The factor could have been a part of Zab{C)- But, unlike for Qab^ we have 

included in Zab only the scaling effects due to interactions. In this case the factor due to 
trivial scaling of the unperturbed operator a has to be added separately. 

Correlation of spins at the cut-off distance a ~ 1 is that of the unperturbed decoupled 
replica models: it is ~ 1 if pairs of spins have same replica indices and otherwise: 

< (7„(0)afe(0)ac(l)arf(l) >~ 5achd + (^aAc (3.8) 

One gets from (3.7): 



< .(OMR) >^ ~ (3.9) 



The matrix ZabU) is assumed to be symmetric in its indices. 



To find < a{0)a{R) it remains to define Zab{0- We shall do it first for the replica 
symmetric fixed point, which is simpler. 

RS. In this case Zab — Z, gab — g, a ^ b. Equation (3.9) takes the form: 

< a{0)a{R) >2 ^ (Z{U)f-^ (3-10) 

Z{^) could be defined from the equations (3.5), (3.6). We shall be deriving correlations 
for the model which is assumed to be already at the fixed point, and not for the crossover 
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behavior when the fixed point is approached. (In any case, for the rephca symmetry 
broken case we have solution for gab for the fixed point only.) In this case we have to 
substitute in cq.(3.6) the fixed point value of g. This has been derived in [1], up to 
second order in e, and reproduced, by a somewhat different technique, in [2]. In the 
normalization that we have chosen at the start of this section, eq.(3.1), the fixed point 
value of g is given by: 

9* = l^ + le' + 0{e') (3.11) 
Prom eqs.(3.5), (3.6), for n = 0, we obtain: 

ZiU)r^exp{^4R} = iRr (3-13) 
We have substituted X = 61og2. Prom (3.10): 

< a{0)a{R) >2 ~ — (3.14) 

with 

= 2A(o) - 7* (3.15) 
By we have denoted the scaling dimension of the operator 

Oab{x) = aa{x)ab{x), a^b (3.16) 

at the replica symmetric fixed point. 

RSB. In the replica symmetry broken case all the matrices are assumed to be in 
Parisi block-diagonal form. One passes to the continuous dependence on matrix indices 
by using the following rules [6]: 

gab ^ git) (3.17) 

{g\a^- tdtg\t) (3.18) 
Jo 

{9\b ^ -2g{t)g - f du{g{t) - g{u)f (3.19) 
Jo 

with 

dtg{t) (3.20) 



14 



and similar expressions for the matrices Z^h and ^oh- The continuous parameter t replaces 
indices, the way described in [6]. It varies originally in the range from 1 to n, natural for 
the matrix indices, but after the limit n ^ is taken the interval of values of t becomes 
[0,1]. The expressions (3.17)- (3.19) are somewhat special since we assume that diagonal 
elements of gah (and of Z^hi ^ah) are zero. Note also that there is no summation aver the 
index a in (3.18), and that for the Parisi matrices all the diagonal elements of the matrix 
are equal, a- independent. 



Using these rules one gets for < (7(0)(7(i?) >^ from (3.9): 



1 " 1 

n- 



= lim^(Z^Ks)). ^ 



The equations (3.3) and (3.4) for Zab and jab take on the following forms: 

^^ = Z{Uh{t) (3.22) 

7(^) = \9{t) + ^{-2g{t)g - du{g{t) - g{u)f) - ^g\t) (3.23) 

We have assumed here that the model is at the fixed point, so that gab and 7^6 are ^ 
independent. g{t) is a fixed point function of t only. It has been found in [3] and it has 
the following form: 

g{t) = ^ ' (3.24) 
[^1, ti<t<l 

^1 = i + + 0{e') (3.25) 

h = 2>gi (3.26) 

We notice that g{t) has a hnearly growing piece for < i < ii, which makes it different 
from the replica symmetric solution: 

9r.s.{t) = const = g^ (3.27) 
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is given by eq.(3.11). We observe that is different in its term from the constant 
part of g(t), gi, eq.(3.25). 

We notice also that the form of g{t) in (3.24) corresponds to the full replica symmetry 
breaking, in the terminology of the spin-glass theory [6]. 

7(t) is defined by the eq.(3.23). It could still be simplified. First, from (3.20), (3.23): 

g^ f dtg{t)^g^-\gl (3.28) 

Keeping the second term of g in the product g{t)g in (3.23) would mean having terms 
~ ~ in this equation, which exceeds the accuracy of our calculations. Second, the 
integral term Jq du{g{t) — g{u))'^ is ~ ~ e^. So it can also be dropped. Then 7(t) takes 
the following form: 

lit) ^{\- j9i)9(t) - ^g'(t) = a^git) - a^g^t) (3.29) 

where 

ai = ^-yilog2~i-^elog2 (3.30) 
a. = 1 (3.31) 

Prom eq.(3.22) 

Z(e,t)~exp{e7(t)} (3.32) 
For the correlation function one obtains from (3.21): 

= J^'dtexp{2^Rj{t)}^^ (3.33) 

Here C,r = logR. 2^R'~f{t) is typically small, < 1. This is the case for example for the 
3-component Potts model and R ~ 10'\ the biggest lattice sizes accessible for numerical 
simulation experiments. In this case it is reasonable to calculate the integral in (3.33) by 
expanding the exponent. In this way one obtains: 

R''^'"' < a{0)a{R) >^ ~ [' dt{l + 2j{t)CR + 2j\t)eR+ ■ ■ ■) 

Jo 

= l + 2^^R + 2^eR + --- (3.34) 
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Here 

7 = / dt'y{t) ~ aig - 025'^ (3.35) 
Jo 

^ = dt^'^U) ~ afg^ - 2aia2^ (3.36) 
Jo 

We remind that the precision of our calculations allows us to keep only the two first 
terms in 7 or in 7^. For g,g^,g^ one finds: 

9^91- I9I (3.37) 
f^gl-2gl (3.38) 
¥^9l-\9t (3.39) 

From (3.35), (3.36), (3.30)-(3.31), (3.37)-(3.39) one obtains: 

1 19 
7 - ^^1 - (log2 + — )c/i 

= ^6-(^log2-i^)6^ + 0(6^) (3.40) 

The expression for the correlation function (3.34) could be given, more conveniently, in 
the form: 

log R^^^°\ a{Q)a{R) >^ = const + 2^^ + - (7)^)^^ + • • • (3.42) 
The term ~ gives correction to the effective scahng dimension of the operator 

Oab{x) = aa{x)ab{x), a^b (3.43) 

A;;^ = 2Af - 7 (3.44) 

-comp. eqs.(3.15), (3.16) for the RS case. 7, eq.(3.40), has to be compared with 7*, 
eq.(3.12). 

The term ~ in (3.42) corresponds to the deviation from the usual scahng form of 
the correlation function. Its coefficient is given by: 

1 27 

2{^'-mc^-g!=-/ + 0{e') (3.45) 
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In numerical simulation experiment it is more appropriate to measure the corre- 
sponding "magnetization", instead of the correlation function. The "magnetization" in 
the present case will be the expectation value of the operator ^ '^a^b Oah = ^ Haj^^b '^aO'b, 
i-G- ^(Sa^b^aft)- This quantity is to be measured for a finite lattice, of size L x L, at 
the critical point of an infinite lattice, and then the finite size scaling analysis is to be 
applied. 

In the theory one obtains < J2a=/=bOab >l, for a finite lattice of size L, from the 
correlation function < J2a^bOab{0) x J2c^dOcd{R) >oc < a{0)a{R) >^ defined on an 
infinite lattice, by putting R — L and taking a square root. In this way one finds, in the 
RSB case: 

\og{L'^'^°'-Q2 OabmL} = log{L'^'°' -J:< aM >L< <7,(0) >a 

a^o a^b 

= consti + + - + • • • (3.46) 

The operator Oab{x) — aa{x)ab{x), in the theory, is a product of local spins for differ- 
ent replicas. In numerical simulations Oab{x) will be a local product of spins for two 
copies of the same random lattice, simulated with different initial conditions for spins. 

Two identical disordered lattices, different starting conditions, they realize in numerical 
experiment two replicas of the theory. We shall have first 

< Oab{x) >L=< aa{x) >l< ab{x) >L (3.47) 

- for a given disorder, and then the average over the disorder, is to be taken. According to 
the replica theory of spin-glasses [6] in the situation when the replica symmetry is broken, 
the space of states is divided into many "valleys" corresponding to different ground states, 
and the summation over replica indices in the thermally averaged quantities corresponds 
to the summation over all these "valleys". Since different initial conditions, in general, 
correspond to different ground states, the summation over various samples (for averaging 
over the disorder) with different initial conditions must correspond to the summation 
over the indices a, b in the quantity (3.47). 
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In the RS case, one obtains from eqs.(3.14), (3.15): 



\og{L'^'^' Y: < O„,(0) >a = log{L2^^°' Yl < >L< a,{0) >r.} 

= const2 + 7*Cl (3.48) 



4 Analyses of the results. Distributions. 

For the purpose of numerical simulation tests we shall give here the numbers for the 
coefficients, for the case of the 3 - component Potts model, e—j^. One obtains: 

3 9 33 
7* ~ -e - (- log 2 - —)e^ ~ 0.050 - 0.019 = 0.031 (4.1) 
8 4 64 

- for the RS case, eqs.(3.12), (3.48). 

7 ~ -e - (^ log 2 - ^)e2 ~ 0, 050 - 0, 024 = 0, 026 (4.2) 
8 4 64 

7^-(7)'^^e' = 0,0005 (4.3) 

- for the RSB case, eqs.(3.40), (3.45), (3.46). 

Characteristic numbers to look at are the values of the products 7^^, (7^ — (7)^)^! 
in eq.(3.46), for L maximal in numerical simulations, L — 10^. One gets, 

7^ = 0.18 (4.4) 

(7^-(7)')Ci = 0.02 (4.5) 

This gives an estimate of ~ 10% on the deviation from scahng, in case of the 3-component 
Potts model. This is only an estimate as we could not know in advance the accuracy of 
the e -expansion calculation. In particular, the scaling violation term (4.5) might still 
be smaller. For this reason we have also done simulations for the 4-component Potts 
model. The results will be presented in the next Section. In the theory, the e - expansion 
values for the coefficients do not make sense for the 4-component model (e = |). In 
particular, the second term in the eq.(3.40) for 7 becomes bigger than the first one, and 
7 becomes negative. On the other hand the model itself should evolve continuously with 
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the number of components q, up to 4 and further. This is because the phase transition 
of the Potts model with random bonds remains second order for g > 4 [9,10,11]. Then 
wc would expect that the effect of the deviation from scaling, if present, should become 
more pronounced as q is increased. So it should be easier observable for the 4-component 
model, compared to the 3-component one. Going from g = 3 to = 4 we go out of the 
perturbative region, where the e expansion is valid, but the qualitative effect of scaling 
violation should increase, if the model is at the RSB fixed point. 

In addition to the above results on the correlation function squared and the associ- 
ated magnetization, one could also look at the corresponding distributions. This implies 
that after the calculation of the thermodynamic expectation values of < aa{x) > and 

< ab{x) >, for two identical lattices, different starting conditions, one changes disor- 
der and performs the thermodynamic measurement again, and gets another values of 

< aa{x) >, < ab{x) > . Having done these measurements many times, one constructs 
the distribution for the values of products of local magnetizations 

< Oab{x) >L=< CTaix) > l< Crb{x) >l (4.6) 

This is instead of summing up the values for products and calculating in this way the 
average over the disorder. [We observe that at the stage of calculating either the distri- 
bution or the average over the disorder of the products (4.6) one could use also the values 
of local products for different points x on the lattice, if they have been measured. But 
we stress again that local products have to be taken first, summation over x second.] 

The distribution obtained in this way, the "overlap function" of local magnetizations, 
could be obtained in the theory from the RG result for the amphtude 

Z(CL,i)~e^W«^ (4.7) 

It is more convenient to study the log of this function: 

Q{t) = log Z{iL, t) = const + ^{t)iL = Qo + l^L (4.8) 

which corresponds to 

Qab = log{(L)2^^°' < aa{x) >L< <Jb{x) > (4.9) 
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We shall define 

q{t) = Qit) - Qo (4.10) 

In analogy with the theory of spin-glasses, to obtain the distribution of values of Qab one 
has to define the inverse function of Q{t), or of q{t), and calculate its derivative. One 
gets: 

q{t) = (4.11) 



, 0<t<ti{^3gi) 



lit) ^ (4.12) 
7i, ti<t<l 

1 1 
7i = am - «2^i = -^91 - (log 2 + —)gl 

= r-(i"'S2-s)^= (4.13) 

We have used equations (3.29)-(3.31) for ^{t) and we have been keeping the accuracy of 
our calculations by dropping extra terms. In particular, in the interval < t < ti we can 
keep linear terms only. This is because the interval itself is ~ e and under integration 
the linear terms become quadratic, ~ e^. One has: 

q{t) = (4.14) 
[qi, ti<t<l 

qi = 7iCl (4.15) 
In the interval < q < qi the inverse function is given by: 

t{q)^^q, 0<q<qi (4.16) 

Its derivative: 

^^Ho<,<„ (4.17) 
dq 

Finally one gets the following distribution function: 

dq [o, g>gi 

Numerical measurement of the distribution N{q) might be complicated. For the 
product of magnetizations on finite lattices of size L, eq.(4.9), the distribution iV(g), will 
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actually be rounded by finite size effects. To have it more distinct one have to look for 
the limit of big L, but then the extra structure in N{q), ior < q < qi, would become 
lower and lower, being of the hight ~ j^. 

The way out could be to look at the distribution of a products of correlation functions 
themselves: 

gi?(i?) = l0g{(i?)^^^°' < aaiO)aaiR) >L< (J6(0)(76(i?) >l} (4.19) 

The index (2) of Q^^jj is meant to tell that we are looking at the overlap function of a 
two-point object. The above considered Qab for magnetizations could have been noted 

In the theory, the corresponding Q'^'^\R,t) shall be given by Z'^{^R,t): 

Q^^\R,t) = log Z\U,t) (4.20) 

= logR, and it is assumed that 1 ^ i? ^ L, L being the size of the lattice. R has to 
be big enough so that the continuum limit theory applies and the fixed point is reached. 
In an analogous way one finds in this case 

dt f ^ + (l-3c/i)5(g-gi), 0<q<qi 
[O, q>qi 

Now, as we increase L, the profile of N^'^\q) will be sharpened while the hight of the 
extra structure ~ remain unchanged. 



5 Simulations. 

In this section, we are going to present some results of numerical simulations that we 
performed in order to check the validity of our results. In particular, we want to try to find 
a way to choose between the two possibles scenarios, replica breaking or replica symmetry. 
The easiest thing that we can simulate, as explained earlier, is the scaling dimension of the 
square magnetization. Thus we have performed the following simulations: On a square 
lattice of size L x L, we simulate two configurations (ex" and a^) of the g-state Potts 
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model with the same disorder, but with different initial conditions. Then we compute 
the product of the magnetization 

Qat = j^ E <<><^-> (5.1) 

Here < a" > means the thermal average of the local magnetization 

= of ■ m« (5.2) 

and m"- is the total magnetization 

^" = 7^ H < (5-3) 
In practice, as we have checked numerically, it turns out that this quantity is the same 



as 



The Hamiltonian of the simulated model is given by 

{iJ} 

where the coupling constant between nearest neighbor spins takes the value 



Jij 



Jq with probability p 
Ji with probability 1 — p 



Measurements were performed on a square lattice with helical boundary conditions. 
Without any lost of generality, we can consider the case where p = \- Then the model is 
self-dual and thus the critical temperature is exactly known. It is given by the solution 



of the equation |13 



1 _ p-PJo 

e-^-^K (5.6) 



1 + (g - l)e-^-^o 

The disorder that we choose to simulate is Jq = 1, Ji = ^. This disorder is in fact quit 
strong in order to avoid problems of cross-over WM. Monte Carlo data were obtained 



by using the well known Wolff cluster algorithm |15|. Due to the strong disorder that 



we considered, we needed to have large statistics over the number of configurations of 
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disorder. Simulations were performed for lattice with size ranging from L — 10 to 
L = 1000. The number of configurations of disorder were 20000 for L = 20-200, 6000 for 
L = 500 and 1000 for L = 1000. For each of these configuration of disorder, measurements 
were taken over ti updates, after to updates for thermalisation. The statistical error SA 
of a quantity A has two contributions, one from the thermal fluctuation, with a variance 
ax, and one from the disorder fluctuation, with a variance a^. Thus the statistical error 
is given by 

where N is the number of configurations of disorder and r is the autocorrelation time. 
For the quantity that we measured, it turns out that the two variance are near equal and 
then 

(Mrc.^(l + ^) (5.8) 

Thus we just need to choose such that ti >> t and we can ignore the thermal fiuctu- 
ations. Then the first step is to measure the autocorrelation time. For the 3-state Potts 
model, we measured r ~ 3 for L=10 up to r ~ 25 for L = 1000. Thus by choosing 
to — ti — 1000, we can safely ignore the thermal fluctuations. For the 4-state Potts 
model, we measured r ~ 5 for L=10 up to r ~ 50 for L — 1000. Again, we choosed the 
parameters to — ti — 1000 except for L — 1000 for which we took = = 2000 and 
this in order to be sure to have thermalized data. 

Results of these measurements for the 3-state Potts model are displayed in Fig. 1. 

In this figure, wc plot ln{L^^'^Qah) versus ln(L). Here, (• ■ ■) means the average over 
the disorder. As explained in section 3, we expected two possible behaviors for this 

quantity, either ln(L^^<^(5a&) — consti -|-7*ln(L) H in the RS scenario (see eq. (3.48)) 

or ln(L^^'^(5a&) — consti + 7ln(L) + (7^ — 7^) In^(-L) + • • • in the RSB scenario (see eq. 
(3.46)). According to the RS scenario, we would expect a scaling behavior ie. a hnear 
(Log-Log) plot, which is in very good agreement with what we obtain. Moreover, we 
estimate that, for the RSB scenario, the deviation from such a linear behavior due to 
the In^(L) term should be of order 10% at L = 1000. We do not see such a deviation. 
Thus our numerical simulation of the 3-state Potts model does clearly favor the replica 
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symmetry solution over the replica symmetry breaking one. 
Performing a fit of the plot 

L^^-^Q^^-L^ (5.9) 

we obtain a value of 7=,, = 0.04 ± 0.002, which is reasonably close to the predicted value 
7, = 0.031 (see eq. (4.1)). 

Fig. 2. corresponds to the same plot but for the 4-state Potts model. As explained in 
section 4, we do not expect that the e-expansion is still valid for g = 4, but we would still 
expect to see a deviation from scaling due to the RSB and we expect that this deviation 
is more pronounced as we increase q. But again, we do not see a deviation from a scaling 
behavior, thus again not seeing any evidence in favor of the replica symmetry breaking 
scenario. Performing a fit of the plot with eq.( ^.9|) , we do obtain here 7=,, = 0.023 ±0.002. 

6 Conclusions and Discussions. 

We consider that the results of numerical simulations presented in the preceding Section 
support the RS fixed point critical behavior. Neither for the 3-component nor for the 
4-component models could we detect the deviation from scaling, characteristic of the 
RSB fixed point. In case of the 3-component model, for which the e-expansion could be 
expected to be reasonably well defined, the value of the slope of the curve for the scaling 
function (3.48) agrees sufficiently well with the value of 7=,, in eq.(4.1). 

Still we would like to make a remark on possible physical significance of RSB fixed 
point, so that it would not appear totally formal. 

The coupling constant Qab in (1-13) is proportional to the overlap function of local 
energies, in a similar way as Zab is proportional to the overlap function of spins. Having 
Qab different for different a, 6, as it is the case for the RSB fixed point, could be attributed 
to a multiplicity of ground states in the statistical model. One starts with different initial 
conditions for spins (one breaks initially equivalence of replicas in this way) and the 
model can end in different ground states. In turn, multiplicity of ground states could be 
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interpreted as a kind of localization phenomenon, in the configurational space of spins. 
Unlike in spin-glasses, in the present model with relatively weak disorder the localization 

would be due to fluctuations, which make disorder important at large distances. It is this 
spontaneous phenomenon that we had in mind behind the notion of RSB in the critical 
model with disorder. 

We would also like to remark that one of the possibilities, why the RSB phenomena 
have not been detected in the particular models considered, could be due to the so called 
"marginal stability" of the RSB fixed point: it is well known that in the linear order 
stabihty analysis of the RSB fixed point (3.24) there exists the so called "zero-mode" 
with the zero eigenvalue. The detailed second order calculations of the stabihty of this 
fixed point shows that to enter the critical regime defined by the RSB fixed point (3.24), 
in general one needs to reach exponentially large spatial scales [16], which could be well 
beyond the lattice sizes of the present numerical tests. 
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A.l The integral of Dp', eq.(2.17). 

^. = /A/<iV^^^3^ (A.1) 

Here Ag = 1 — |e. We change the variable y' : y' — yt. This gives: 



2/3 



i<\y\<a 



.1.6. r(i + a)r(i + /3)r(-i-a-/3) 
^ -'^^ r(-«)r(-/3)r(2 + « + ^"^-'^ 

Here a — — ^ ~ ~| + f ^' ~ = — 1 + |e, and we have used the known result: 

i ' ' ' ' r(-a)r(-/5)r(2 + « + ^ ' 

For its derivation see e.g. [12]. Putting the values of a, (3 and expanding in e one obtains: 

= ^(l + e61og2) (A.4) 



Finally one gets: 



ii: = 61og2. 



I, = 2^{l + eK)^{o'^-l) (A.5) 



A.2 The integral of Df\ eq.(2.19). 

j d^y J d^y' < a (0) e {y)e{y') a (oo) >< e{y)e{y') > (A.6) 

This integral has already appeared in the calculations in the paper [2]. It is finite, the 
limit of e 0. This could also be seen without calculations, by using the operator 
algebra. In fact, when y — > one has: 

<^)<y) = ^^^0) + • • • (^-^^ 

As Ag 1, the integration over y around is finite. 
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When y, y' — > 0, y' ^ y, one has: 

a{0)e{y)eiy') = , .^^^ ^ a(O) + • • • (A.8) 

\y\ \y I 

Again, the integration around is finite. 

Finally, for the configuration of — > y one obtains: 

e{y)e{y') = ^(1 + ■ ■ ■ + \y - y'\'T{y)f{y) + •••) + ••• (A.9) 

T{y),T{y) are components of the energy- momentum operator. In the OPE (A.9) the 
term ~ TT is the first after 1 which contributes when integrated in (A. 6). Some extra 
terms, which vanish under integration, like {y — y')'^T[y) or [y — y')'^T{y) have been 
dropped. Putting fti 1, one gets 

< a{0)e{y)e{y')a{^) >< e{y)e{y') > (A.IO) 

=< a(0)a(oo) > ^^^(1 + \y- y'\^ < <j{0)T{y)f{y)a{oo) >+...) + ... 

The first term, ~ l/\y — y'\^, leads to quadratic divergence in the integral (A. 6). It is 
automatically subtracted in the analytic e-expansion calculation. It corresponds, in fact, 
to the constant shift renormalization of the action of the theory which is irrelevant. The 
next term is finite, when y' — > y. So the integral (A. 6) will be finite. 

A.3 The integral of Df \ eq.(2.21). 

73 = / A' / A(< a{Q)e{y)e{y')a{:^) >f (A.ll) 

The variable y' could be scaled out by using invariance of correlation functions w.r.t. 
global dilatation. In general, when a set of operators Oi, O2, • • • , 0„ is projected on the 
operator On+i, placed at infinity, one has: 

< Oi{xi)02{x2) ■ ■ ■ On{Xn)On+l{oo) > 

= A^i+^-"+^"-^"+i < Oi(Axi)02(Ax2) • • • 0„(Ax„)0„+i(oo) > (A.12) 
where A is a dilatation parameter. For < a{0)s{y)e{y')a{oo) > this gives: 

< a{0)6{y)e{y')a{oo) >= A^^^' < a{0)6{\y)6{\y')a{oo) > (A.13) 
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Putting \ — l/y', one obtains: 

< a{Q)e{y)e{y')a{^) >= \y'\-^^^ < a(0)£(^)£(l)a(oo) > (A.14) 

Using this relation in the integral (A. 11) and changing the variable y, y — y'y, one 
obtains: 

73 = / d'y'lyf-'^^ J d'y{< a{0)s{y)e{l)a{oo) >)' (A.15) 
Putting = 1 — |e and integration y' between the cut-offs 1, a, one gets: 



73=/ £y'\y\-'+'^h = 27rha'' - l)h (A.16) 

Jl<\y'\<a be 

h = j (fy{< a(0)£(y)£(l)a(oo) >f (A.17) 

In the Coulomb gas representation the correlation function < (reea > could be presented 
in the following form: 

< (j(0)£(t/)£(l)a(oo) > 

J d2ii|y|^"'^"^|y-l|^"'|ii|^"'''*+|y-ii|^"^"+|l-ii|^w (A.18) 
Here is the normalization constant. It can be fixed by using the operator algebra: 

r ^ 0, a{0)a{r) = + ■ ■ ■ (A.19) 

r ^ 0, V^MVaAr) = -^V2aM + ■ ■ ■ (A.20) 
For the correlation function 

< a{0)a{r)s{y)e{y') > 

^Njd'uK v^MVaAr)VaMVaAy')VaM > (A-21) 

This gives, in the limit of r ^ 0: 

< e{y)e{y') > = N-^ f d\ < F2ao(0)K.(l/)K.(l/')K+(w) > (A.22) 



1 2 A,. 

1 



X y d'ii|ti|^w \u - yf'^''^^\u - w (A.23) 
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Changing the variable u in the integral: 



y-y 

one obtains, after simple algebra: 

1 



u = X ^ (A.24) 



N- 



]-^l£u\ur<\u-ir< (A.25) 



Here 2A^ = 4Ai2 — 4Qii2 — 8q;i2Q;o- Using = —a+/2, 2ao = a+ + a-, a+a- = 1, one 
gets also 2A^ = -2 + Sal- Prom (A.25): 



N 



Here = | — e. One obtains: 



(y'd\iMr2"+iM-ir2"+)-^ (a.26) 



_ r^K)r(2-2c.^) _ r^(|-e)r(-| + 2e) 

P(l-a^)r(-l + 2a^)-^ P(_l + e)r(|-2e) ^""-'^^ 

We have used again the integral (A. 3). In our further calculations the normalization 
constant N will multiply expressions with singularities ~ 1/e, but not This implies 
that we don't need to keep the e correction of N. For e = 0, and using the usual properties 
of F-functions, one finds from (A. 27): 



N = 



2F3(-|) _ 2 T\-l: 



97rF3(-i) 9sin(7r(-|))r(|)F3(-i) 

2 r^(-i) 



(A.28) 



or 



Returning to the integral Is in (A. 17) and using the expression (A. 18) for the correlation 
function < aeea >, one gets: 

is = N'^I (A.30) 

7 = Jd^yJ d\ij d\2\y\^°'"''\y - l\^'''\ui\^"''"+\ui - 1|^w|mi -yj^w 

x\u2\^^^''+\u2 - 11''"^"+ 11*2 - (A.31) 
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The integral I is calculated in the Appendix B, with the result: 

71 



I^-^7' + 0{e) (A.32) 



7 is defined in (A. 29). For Is one gets: 



73 = A^'/ = -^ + 0(e) (A.33) 
and then for the integral I^, (A. 11), one obtains from (A. 16): 

+ (A.34) 

6 be 

This expression is in fact incomplete. It has to be corrected. There is an extra term 
in it, ~ 1/e^, which is missed in the analytic technique of calculating the integral I in 
the Appendix B. It could be recovered by using the original expression for the integral 
Is in the eq.(A.ll) and the operator algebra. Configuration which is responsible for the 
leading singularity in is either 

y,y'^o (A.35) 

or 

i<lz/'l<|y|, y,y'^o (AM) 



For (A.35) one has, by operator algebra: 



a{0)e{y)e{y') = ^aiO)e{y') + 

\y\ ' 



\y\ ' \y 

D — D^^ is the OA coefficient. For the integral /a, this gives: 



jcPy{< a{Q)e{y)e{y')a{^) >f 

= / 27^\y'\d\y'\ j 27r|y|%|--^^^ + 

Ji<\y'\<a J^<\y\<\y'\ \y\ ^\y I 

^Att'^D^ [ d\y'\\y'\-'+^' [ d\y\\y 

Jl<\y'\<a il<|2/|<|2/'| 



l+3e _^ 



An'D' [ d\yr'^'^^\yT + --- 

Jl<\y'\<a 3e 

27r'L>2-i-a^^ + • • • (A.38) 
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We have used: < a{0)a{oo) >— 1, 2A^ = 2 — 3e, and we have kept the leading term 
when integrating: ly']-^^ — 1 ^ \y'\'^'^, a^^ — 1 ~ a^^. 

Adding to (A. 34) the piece in (A. 36), multiphed by two because of two equivalent 
configurations, (A. 35) and (A. 36), one finally obtains: 

^3=(47r^^'^-^)a^^ + 0(l) (A.39) 

The leading singularity, ~ is missed in the calculation of the integral / in the 
Appendix B for the following reason. When 73 is put in the form of (A. 16), (A. 17), the 
factor 1/e is already in front, 1/e^ singularity of Is would be produced by 1/e singularity 
of /a, or /, eq.(A.31). Instead we find the result (A. 32) for /, / ~ 1. In the integral 
/3,(A.17), or in /, (A. 31), the integration over y is performed over the whole infinite 
plane. This is correct to define the finite piece of the integral, but the singularity 1/e get 
cancelled in this way. 

In fact, let us consider the integral I^, (A. 17). There are two configurations which 
lead to 1/e singularity in the integral (A. 17): 

< |y| < 1 (A.40) 

and 

K ll/l < 00 (A.41) 

For the first one we use: 

'^(^Hv) = 7:^^^(o) + • • • (A-42) 



\y 



and we get: 



< a(0)£(y)£(l)a(oo) >= < cr(0)£(l)(7(oo) >= (A.43) 

\y\ ^ \y\ ^ 



We have used < a{0)e{l)a{oo) >— D. For the second configuration, (A.41), we use: 

a{0)e{l) = Da{0) + ■■■ (A.44) 

and we get: 

< a{0)6{y)e{l)a{oo) > = < a{0)e{l)e{y)a{oo) > 
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^ D< a{0)e{y)a{oo) > 

= < ^(OMlMoo) ^ (A.45) 

In the last line we have used the scaling properties of the function < a{0)e{y)a{oo) > . 
In the integral I3, (A. 17), the above two configurations will provide the following contri- 
bution: 




(A.46) 



Next we change the variable in the second integral: \y\ —>■ l/\y\. Then we obtain: 

/3 = 27r/ d\y\\y\-'+'' + 27i [ d\y\\y\''-^' + ■ ■ ■ (A.47) 

JO<|j/|<l J0<|2/|<1 

If the first integral is ~ 27r/3e (by extending the integration up to 1), then, by analytic 
continuation in e, the second integral should be defined to be ~ — 27r/3e. In the result 
the two contributions ~ 1/e cancel one another. This is what is happening with the 1/e 
terms in the calculation of the integral / in the Appendix B. The calculation there is all 
based on the analytic continuation technique. 



B Appendix 

The integral /, eq.(A.31). The integral is of the following general form: 

I ^ jdH Jd'x Jd'y\t\'^'\t-l\'^' 

x|x|2"|a; - l\'^'^\x - t\'^''\y\'^''\y - l\'^^\y - (B.l) 

We shall calculate it for the values of exponents a', /?', a, p which will be specified 
below. By the standard technique, see e.g. [12], this integral could be factorized into the 
following sum of products of contour integrals: 

I = -{A''\s((3')s\(3)ji-^ + s{(3')s{(3)s((3 + p)^-^ + s{(3')s\(3 + p)j'f^\ 
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+jt\siP' + 2p)s\P)A-^ + + + siP')s\P)A-^]} (B.2) 

Here = sinvr/^, etc. The contour integrals are defined in the following way: 

= ('dtfdxfdy{tr\l-tf 
Jo Jo Jo 

x{xni - xf{t - xYivni - yf{t - vY (B.3) 

= 2 /'dt j'dx fdy{---) (B.4) 
JO Jo 

= f'dt f'dx f'dyi---) (B.5) 
JO Jt Jt 

The symbol (■ ■ ■) in (B.4), (B.5) stands for the same expression as in (B.3) except that 
the variables are put in the order corresponding to the order of integration, so that the 
differences of variables are always positive. E.g. the factors {t — xY, {t — yY in (B.3) will 
be in the form {x — tY, {x — tY in the integral (B.5). 

/OO ft ft 

dt dx dy{tr'{t-lf 

X {xYix - lYit - xYivYiy - ^)^{t - vY (b.6) 

"^~> ^ 2J^ dt dx dy{- ■ •) (B.7) 



/oo poo poo 

dt dx dy{- ■ •) (B.8) 



The integrals j| can be put in the same form as the integrals j} ,^2 ,js 

under the change of variables t ^ 1/t, x ^ 1/x, y ^ 1/y. One obtains: 

A"^ = C dt dx j' dy{tf\l-tf 
Jo Jt Jt 

X {xfii - xf{x - tYiyfii - yf{y - tY (B.9) 

A"^ = 2 f'dt f dx C dy{---) (B.IO) 

Jo Jo Jt 

■■ dt t dx t dy{- ■ ■) (B.ll) 
Jo Jo Jo 
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Here 

a' = -2 - a' - - 2p, a^-2-a-(5-p (B.12) 
Originally the values of the exponents are given by: 

a' = 4a,a, = i-^ (B.13) 

(3' = Aal^^ + e (B.14) 

a = 2aaa+ = + | (B.15) 
4 

/3 = 2q;^q;+ = -- - e (B.16) 
4 

p = 2aea+^---e (B.17) 

Here cKo-, correspond to the Coulomb gas operators 

^aaix) -. exp{iQ;^(/7(a;)} : ; Va,{x) -■ eicp{iaeip{x)} : 

which represent the operators of spin a{x) and energy e{x)] a+ corresponds to the screen- 
ing operator. In particular: 

a, = ai,2 = (B.18) 
For the spin operator, if is put in the form: 

(which corresponds to the Kac table of the size (2p — 1) x (2p — 2)) then the Potts model 
spin operator is 

'^-yp,P-i,V^ (B.20) 
Vp^p-i ix) is the conjugate Coulomb gas operator (with respect to Vp^p-i{x)), which we 
are actually using: 

V^(^) exp{ia^<^(x)} : (B.21) 

1 + p p 

a- H — ( 

2 2 

One gets the values of exponents in (B.13)-(B.17), using in addition the relation between 



'^p^iTT = = — - — tt- + -a+ (B.22) 



p and 



Pi< - 1) = ^ (B-23) 
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which follows from (B.19), and the usual relation q;+q;_ = —1 for the Coulomb gas 
parameters. 

We have found that the calculation of the integral / is simpler for the values of 
exponents a', a, /3, p obtained after the transformation of variables: 

t^l-- x^l--, y^l-- (B.24) 
t X y 

One gets then the following values: 

a' = l + e, f3' = -l + ^e (B.25) 

a^~-e, /3 = l + ^e (B.26) 
p=-^-e (B.27) 

The exponents a', a of the contour integrals J2~\ Js '* ^qs. (B.9)-(B.ll), are obtained 
from (B.12): 

a' = i-| (B.28) 

a = -i + I (B.29) 

We shall give next some details on the calculation of the contour integrals , j2^^ , , 
which enter into the decomposition of the integral /, eq.(B.2). The calculation of the 
integrals j3~\j2~\ji~^ is respectively analogous. 

= dt - tf t dx - xf{t - xY 
Jo Jo 

X fdyy^{l-y)^{t-yy (B.30) 
Jo 

We change variables: x — xt, y — yt and next we drop the tildes of the new variables 
X, y. 

Jo Jo 

['dyy^{l-yy{l-ytf (B.31) 
Jo 
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X 



We expand next the factors (1 — xt)^, (1 — yt)^ : 



^i_^tr^f:^-^{xtr (B.32) 



fe=0 



Here = (-/?)(-/? + 1) ■ ■ ■ {-f3 + k - 1). One obtains: 

hi k2\ Jo 



X dx x'^+'^Ui - xY C dy y'^^'^Hi - vY (b.ss) 

Jo Jo 

Next we use the integral 

Jo ^ ^ r{2 + a + b) ^ ^ 



to obtain: 



•(+) ^ (-/g)fc2 r(3 + a' + 2a + 2p + ki + k2)T{l + p') 

^i' r{4 + a' + P' + 2a + 2p + ki + k2) 

r{i + a + k^)r{i + p) r(i + g + k^W + p) 

r{2 + a + p + ki) r{2 + a + p + k2) ^ ' ' 

By using repeatedly the recurrence relation for the F-function, Y{z^ 1) = zY{z\ one can 
put (B.35) into the following form: 

jt^ = it'^st' (B.36) 

(+) ^ r(3 + + 2^ + 2p)r(i + /jQ r(i + a)r(i + p) 

r(4 + a' + + 2a + 2p) ^ r(2 + a + p) ^ ^ ' 

^(+) ^ ^ ^ (-/^)fci (-/^)fc2 (3 + + 2a + 2p)fc^+fc, 



iti 2 



/ci! /cs! (4 + a' + + 2a + 2p)fe,+fe2 

^^1±^J1±^ (B.38) 
(2 + a + p)fei (2 + a + p)^. 

Substituting the values of exponents (B.25)-(B.27) one obtains: 

- r(-i-|.) ( r(-i-2e) ' 

c.(+) ^ y^ y^ ~ ~ i^)fe2 (-| - ^)fci (-| - e)fc2 (-1 - 3e)fci+fc2 {jx, Aa\ 

' XX (-1 - 26k (-1 - 26k (-1 - fek^,, ^ ■ ^ 
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We shall do calculations by expanding in e. For the renormalization group equation we 
need to keep the first two terms only. 

For 7i+) in (B.39) a simple calculation gives: 

^'^^-te^^-l'-^"') + 0{e) (B.41) 



We have defined here: 



in-w 

r(-i) 



7 = (B-42) 



-^(4)-^(-^)-^ + ^ (B.43) 
il){z) is the standard ■0-function: 

^{z) ^ j-^\ogT{z) (B.44) 

The calculation of S^'^ is more involved. First we check the convergence of the series 
in (B.40). In general one has: 

^p.(fcr^ (B.45) 

Using this asymptotic form one gets the following estimate for large ki, k2 in (B.40): 

^i'*"^-^EE(^i)~^~*(M~^~*(^i + M"^' (B.46) 

fei k2 

(We observe that A;! = (l)fc)- So the series converge and we can do safely its e-expansion. 

Using the specific values of the parameters in the series (B.40) one can develop it in 
the following way: 

oo oo 

(0,0) + 2(0,1) + (1,1) + 2 Y: (•••)+ 2 E (•••) (B-47) 

fc2=2(fci=0) k2=2(ki=l) 

(0,0) stand for the term of the series ki = 0, /c2 = 0, etc. We have dropped the part of 
the series corresponding to (/ci = 2, • • • oo, A;2 = 2, • • • , oo) since it is ~ e^. We are keeping 
only terms ~ 1 and ~ e in S['^\ Explicitly one gets: 

~ ' + ' 1 (-i-2e)(-l-i.) ^^-''^ 
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1 ^-|-2e^-l-i6)(-|6) 

^ (A; + 2)! (-1 - 26W (-1 - i6W2 

oo 

+2E 



-1 - |e) (-1 - le)k+2 (-1 - e) (-| - e)k+2 (-1 - 36)fc+3 



,=0 1 (^ + 2)! (-1 - 2e) (-1 - 2e),+2 (-1 - feWs 

To simplify further, we use the following relations: 

{a)k+2 = a{a + l)(a + 2)^, 

{a)k+3 = a(a + 1) (a + 2) (a + 3)fe (B.49) 

In particular: 

(-1 - le),^2 = (-1 - - ^e), ^ ^6(1), = ^ek\ (B.50) 

(-^ - e)k+2 ~ (-^).+2 = (B-51) 
etc. Using these simple rules, after some algebra one gets Si^^ in the following form: 

S[-'^ = l + e{3s-^-) + 0{e') (B.52) 

where 

It remains to calculate this sum. We shall use the following general results: 

{a)k{b)k ^ r{c)r{c-a-b) 

to k\ic), ric-a)r{c-b) ^ ■ > 

- k\ ib),_c-l c-1 c-b 

to{k + 2)\{c),-b-l^c-2^b-2''^ ^^-^^^ 

Here 

/c = ^(c-6) -^(c-2) (B.56) 

Eq. (B.55) could be derived by using the sum (B.54), which is standard, and some rather 
simple algebra. 

From (B.55) and (B.53) one gets in a straightforward way: 

« = -| + f (B.57) 
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The constant k is defined in (B.43). By eq.(B.52) 



« ^ + e(3. - ^) = 1 - + ^e«: (B.58) 

Returning still back to eqs. (B.36), (B.41) one obtains: 

= 7i^^^i+^ = 1^(1 -\e- Ke) + 0{e) (B.59) 

The integral eq.(B.ll), is of the same form as with only the exponents a', a 
replaced by a', a, eqs. (B. 28), (B.29). The calculation follows the same lines, with the 
result: 

jt^=f(-| + f) + 0(6) (B.60) 

,{+) A-) 

For calculation of ^2^^ it is useful to use the following linear relation of the integrals 

U2 = — , ^ . ui ——^U2 (B.61 

Here we have redefined, for the purpose of this particular calculation only, 

«2 = i#\ «i=jf^ (B.62) 

and U2 is a new integral: 

f\ pt poo 

U2= dte\l-tf' dyy''{l-yf{t-y) dx x'' {x - if {x - t)" (B.63) 
Jo Jo Jl 

s{a) in (B.61) is sinTro;. The relation (B.61) is obtained in a standard way by doing 
transformation of the contours of integration [12]. 

Putting the problem this way, to calculate j^^^ we have to calculate the integral M2, 
(B.63). Because of specific values of the exponents the calculation of the integral U2 is 
simpler, as compared to the direct calculation of jg^^ 

We do first the change of variables in (B.63): x —>■ 1/x, y — > ty. One gets: 

U2 = [' dt t'+-'+-+P{l - tf f' dy |/"(1 - yy{l - tyf 
Jo Jo 

X dx x^(l-xf(l-xtY (B.64) 
Jo 
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a was defined in (B.12), (B.29). Next we expand the factors (1 — ty)^, (1 — xty and 
proceed like we did for the integral starting with eq.(B.31). We get U2 in the 
following form: 

U2 = I2S2 (B.65) 

_ _ r{2 + a' + a + p)r(i + f3') Tji + c.)r(i + p) r(i + 5)r(i + 

r{3 + a' + p' + + p) r(2 + a + p) T{2 + a + p) ^ ■ ^ 

g ^sp^ i-f^)k{-p)i + {2 + a' + a + p)k+i 

' Y z ^' ^' {2 + a + p)k{2 + a + P)i{3 + a' + f3' + a + p)k+i ^ ' ' 
Substituting the values of the parameters, eqs.(B.25)-(B.29), one obtains: 

_ r(i - e)r(|e) r^H - e) r(| + |)r(2 + |e) 

r(| + |) r(-i-2e) r(| + 2e) ^^-^'^ 

' rr ^! ^! (-|-2e),(| + 2eM| + tW, ^^"''^ 
For 72 one gets: 

72 = ^7(1 + (f - - 2e« - 3eV'(-^) + SeV'll)) + 0(e) (B.70) 
For S2 one gets asymptotically, for large k^l: 

- the series is convergent. In order to develop in e, S2 can be decomposed in the following 
way: 

-52 = (0,0 + (1,0 + (^ + 2,0 (B.72) 



or explicitly: 

-§2 ~ E 



1=0 ^' (I + 2e)/ (1 + f 



^ 1 ' (-i-2^)(i + i)S ^! (l + 26MI + l)^ 

2^\-|)(|),^„(^ + 2)!(|),^ U (I). ^^-'^^ 

In the last term we have put e = because of the factor (— |e) in front. Finally one gets: 

13 3 
S2^Si--{l--e)s2 + -essS4 (B.74) 
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with 

°° (-) (-) 

By using the results for sums in (B.54), (B.55) the calculation of Si, S3, S4 is straightfor- 
ward. One obtains: 

5 6 27 1 

si = --(l + (--y)e + «:6 + 36V^(--)-3e^(l)) + 0(e2) (B.79) 

«3 = -\ + ^ (B.80) 

S4 = -^0' (B.81) 

For the sum S2, we observe at first that 



I 



(B.82) 



So one gets: 



_ .2 6 - (1 + 6), (1-6), 1 

- ^3 + 2)g^(|T^iTiT] ^^-^'^ 



i + f ^(l + e), (|-e), 1-6 + Z 



— e 



E 



2 ^\ Zl^ 71 I o^N 2 



/! (| + 26),| + f + / 



e^; /! (i + 2e),^ 2 | + f + /' 



3 ^ /=o 

Calculation of the first sum is straightforward. One obtains: 
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The second sum could be calculated by decomposing 

1 111 



V3 ' VV3 

and using the formulas: 



(B.85) 



^(a), 1 _ r(i - a)r(b) 
One gets in this way: 

From eqs.(B.83), (B.84), (B.88) one obtains: 

Putting together the results for si, S2, S3, S4 into eq.(B.74) for 5'2 one gets: 

^2 = -^7(1 + + + SeV'l-^) - ^ei^il)) + 0(e2) (B.90) 

Returning back to the eqs.(B.65), (B.68) one derives: 

{is = -^(1 -7e + 2Ke) +0(e) (B.91) 

Going still back to the relation (B.61), putting the values of the parameters and substi- 
tuting the value of ui obtained earlier, one finally derives: 

) = 2u2 = ^(^ - 3/^) + 0(e) (B.92) 

The integral j2~\ eq.(B.lO), is of the same form as js^"* with only the exponents a', a 
replaced by a', a, eqs.(B.28), (B.29). The calculation follows the same line. There is one 
additional detail in the calculation. In the decomposition of the sum 5*2, corresponding 
to the integral j^"*, one gets, in particular, the sums: 

_ v^ (| + g)K-i-e)/(i-e); .pjqoA 

*l-2^ 71 r5 , 6W2 , e\ l^-y-^J 
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The sum si is transformed as follows: 
and next relation is used: 

(fe^)^ = + 0{^) (B.96) 

So, up to linear order in e, one obtains s\ in the form: 

In the rest of the calculation one proceeds in the way analogous to the calculation of the 
sum 

S2 for the integral j2^\ eq.(B.83). The sum S2, of the integral j2~\ eq.(B.94), is dealt 
with in a similar manner. Finally one obtains the following result for j2~^: 

A-^^-^^{l-le-Ke)+0{e) (B.98) 

■73 )Jl ■ 

We redefine: j^^^ = u^, and we shall use the following linear relation between the 
integrals: 

Us = A\i + 2ABu2 + (B.99) 

Here Ui = U2 is defined in (B.63), u\ is a new integral: 

fi 



dtr{l-tf dx x^{x-l)l^{x-tY dyy"{y-lf{y-t)P (B.lOO) 



A and B in (B.99) are coefficients: 

A = B = -'(^ + f^ + P) (B.lOl) 

s{a + p) s{a + p) 

The relation (B.99) is obtained by transformations of the contours of integration [12], 
similar to the relation (B.61) 
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The integrals ui, U2 had already been calculated. To define 1*3(^3^^) we have to define 
ul, which turns out to be easier. 

We change the variables: x 1/x, y 1/y. This gives: 

ul = C dt - tf C dx x^il - xfil - xtY C dy y^(l - yf(l - ytV (B.102) 
Jo Jo Jo 

a = —2 — a — (3 — p. Next we expand the factors (1 — xty, (1 — yty. This leads to the 
following form of u\ : 

ul = 7i5i (B.103) 

, ^ V{l + a')V{l + (3') V{l + a)T{l + (5) ^ 
r(2 + «' + /?') ^ T{2 + & + 0) ' 

_ r(| + 6)r(|6) r(| + |)r(2 + |6) 

- r(| + ie) ^ r(| + 2e) ^ ^^-'""^ 

S* = ^^ (-P)ki-P)i (l + «)fc (1 + Q:); jl + a')k+i 
' Y i ^' ^- i2 + a + f3)ki2 + a + f3)ii2 + a' + P%+i 

_ (I + (I + (j + (1 + D' (3 + min"^^ 

- 4^4^ ^! /! (| + 2e).(| + 2eM| + i6).H-, ^^-'"'^ 

One easily checks convergence of this sum. 

It turns out that the contour integral ^3^"* , like also j[ ^ , enters into the decomposition 
of the integral /, eq. (B.2) , with a coefficients ~ e^. This immediately follows from eq. (B.2) 

after the substitution of the values of the exponents /9' = — 1 + |e, /9 = 1 + |e p = — | — e 
eqs.(B.25)-(B.27). As a result, it is sufficient to define the leading ~ 1/e parts of the 
integrals 



One finds fi-om (B.104): 



7i* = |(^f + 0(1) (B.106) 



and from (B.105): 



This gives: 



St - (E||^)= = (^)^ - (5)V + 0(.) (B.107) 



'''i = f'iSl^^ + 0{l) (B.108) 

be 
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Finally from (B.99), one gets: 
or 

= 0(1) (B.llO) 
This implies that ^3^'' will not contribute to / (to the finite part of /, to be more precise). 
Similar calculation for the integral gives: 

jf^ = -^ + 0(l) (B.lll) 

We come back now to the eq.(B.2). Substituting the values of the exponents in 
the coefficients, developing the coefficients in e and keeping the leading terms only, one 
obtains: 

+jt\^'s'^A-' - ^'e'^A-' - vrVji-)]} (B.112) 
Here we have noted ^e — e. The leading terms of the contour integrals are the following: 

jf ^ ^ iJ, jf ^ ^ y (1 - 3/.), jf ) ^ 1 (B.113) 

A-^ - A'^ - A-^ - + f ) (B.114) 

Taking them into account, the expression (B.112) for / can further be reduced: 

It turns out finally that we need to know only the leading terms of the integrals. Still, 
the leading terms for A'^^ ^^'^ A~^ finite terms. So we had to calculate them, the 
singular ~ ^ terms for all the contour integrals would not be sufficient. 

Substituting the values (B.113), (B.114) of the integrals, replacing back e by |e and 
using K— + |, one finds after some simple algebra: 

7 = -^7' + 0(e) (B.116) 
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